LESSONS FROM B Xg^i IN TWO HIGGS DOUBLET MODELS 

F. M. BORZUMATI 

Laboratoire de Physique Mathematique et Theorique, Universite Montpellier II, 
F-3405 Montpellier Cedex 5, France 

C. GREUB 

Institut fiir theoretische Physik, Universitdt Bern, 
Sidlerstrasse 5, 3012 Bern, Switzerland 



ABSTRACT 

The next-to-leading order predictions for the branching ratio BR(_B Xs^) are 
given in a generahzed class of two Higgs doublets models. Included are the recently 
calculated leading QED corrections. It is shown that the high accuracy of the Stan- 
dard Model calculation is in general not shared by these models. Updated lower limits 
on the mass of the charged Higgs boson in Two Higgs Doublet Models of Type II are 
presented. 



Talk presented by C. Greub at the 
XXIX International Conference on High Energy Physics 
Vancouver, B.C., Canada, 23-29 July 1998 



LESSONS FROM B Xg^i IN TWO HIGGS DOUBLET MODELS 



F.M. BORZUMATI 

Laboratoire de Physique Mathematique et Theorique, Universite Montpellier II, F-3405 Montpellier Cedex 5, France 

E-mail: francesc@lpm.univ-montp2.fr 

C. GREUB 

Institut fur theoretische Physik, Universitdt Bern, Sidlerstrasse 5, 3012 Bern, Switzerland 

E-mail: greub@itp.unibe.ch 



The ncxt-to-leading order predictions for the branching ratio BR{_B — > Xg'y) are given in a generalized class of two Higgs doublets 
models. Included are the recently calculated leading QED corrections. It is shown that the high accuracy of the Standard Model 
calculation is in general not shared by these models. Updated lower limits on the mass of the charged Higgs boson in Two Higgs 
Doublet Models of Type II are presented. 



1 Introduction 

Two Higgs Doublet Models (2HDMs) are conceptually 
among the simplest extensions of the Standard Model 
(SM). They contain additional sources of flavour change 
due to their extended Higgs sectors. Studies of the 
B Xsj decay in this class of models, therefore, can 
already test how unique is the accuracy of the SM re- 
sult for the branching ratio BIU_B — > Xgj) at the next- 
to-leading (NLO) order in QCDI^I, even requiring calcula- 
tions at the same level of complexity as the SM one. They 
can obviously provide also important indirect bounds 
on the new parameters contained in these models. In 
spite of their apparent simplicity, indeed, 2HDMs have 
not been-cprrectly constrained in ongoing experimental 
searches 013. 

The well-known Type I and Type II models are par- 
ticular examples of 2HDMs, in which the same or the 
two different Higgs fields couple to up- and down-type 
quarks. The second one is especially important since 
it has the same couplings of the charged Higgs H'^ to 
fermions that are present in the Minimal Supersymmetric 
Standard Model (MSSM). The couphngs of the neutral 
Higgs to fermions, howev£i_ have important differences 
from those of the MSSMBI. Since, beside the W, only 
charged Higgs bosons mediate the decay B ^ Xs'y when 
additional Higgs doublets are present, the predictions of 
BR(B Xs"f) in a 2HDM of Type II give, at times, a 
good approximation of the valucLof this branching ratio 
in some supersymmetric models □. 

It is implicit in our previous statements that we 
do not consider scenarios with tree-level flavor chang- 
ing couplings to neutral Higgs bosons. We do, however, 
generalize our class of models to accommodate Multi- 
Higgs Doublet models, provided only one charged Higgs 
boson remains light enough to be relevant for the pro- 



cess B Xsj . This generalization allows a simulta- 
neous study of different models, including Type I and 
Type II, by a continuous variation of the (generally com- 
plex) charged Higgs couplings to fermions. It allows also 
a more complete investigation of the question whether 
the measurement of BR(i3 Xgj) closes the possibility 
of a relatively light not embedded in a supersymmet- 
ric model. 

This summer (1998), a new (preliminary) measure- 
ment of this decay rate was reported by the CLEG Col- 
laboration! BR(B Xs-f) = (3.15±0.35±0.32±0.26) x 
10~^ , which, compared to older results, is based on 53% 
more data (3.3 x 10® BB events). The upper limit al- 
lowed by this measurement, reported in the same paper, 
is 4.5 X 10-4 at 95% C.L.. The ALEPH Collaboration has 
measuredB BR(B ^ Xs^f) = (3.11 ± 0.80 ± 0.72) x lO"" 
from a sample of b hadrons produced at the Z-resonance. 

The theoretical prediction for BR{B Xg^) has a 
rather satisfactory level of accuracy in the SM. The main 
uncertainty, slightly below ±10%, comes from the ex- 
perimental error on the input parameters. The more 
genuinely theoretical uncertainty, due to the unknown 
value of the renormalization scale Hb and the matching 
scale fiw, which was unacceptably large at the leading- 
order (LO) level, was reduced to roughly ±4% when 
the NLO QCD corrections to-the partonic decay width 
r(6 Xg"/) were completed 13. (See ref. El for the mile- 
stone papers which brought to the complete LO calcu- 
lation.) In addition, non-perturiiative contributions to 
BR{B — > Xsj) , scaling like 1 /m^H and 1 / m^EI, were com- 
puted. Very recently, the leading QED and somfijdaases 
of electroweak corrections were also calculatedli^'£i O. 

Following the procedure describod in refs ]M3, and 
including QED coirections as in rcf. t3, we obtain a theo- 
retical predictionEj in agreement with the existing data: 

BR(B Xs-f) = (3.32 ± 0.14 ± 0.26) x 10"^. (1) 



The first error in (|^) is due to the fib and scale un- 
certainties; the second, comes from the experimental un- 
certainty in the input parameters. 

A detailed study of B ^ X^j at the NLO in QCdS 
in 2HDMs, on the contrary, shows that the NLO correc- 
tions and scale dependences in the Higgs contributions 
to BR(i? — ^ Xg'j) are very large, irrespectively of the 
value of the charged Higgs couplings to fermions. This 
feature remains undetected in Type II models, where the 
SM contribution to BR{B Xgj) is always larger than, 
and in phase with, the Higgs contributions. In this case, 
a comparison between theoretical and experimental re- 
sults for BR(_B Xs'^) allows to conclude that values 
of m^± = 0{mw) can be excluded. Such values are, 
however, still allowed in other 2HD models. 

These issues are illustrated in Sec. 4, after defining 
in Sec. 2 the class of 2HDMs considered, and presenting 
the NLO corrections at the amplitude level in Sec. 3. A 
brief discussion on the existing lower bounds on m^± 
coming from direct searches at LEP is included in SeG,[3, 
This rests on observations brought forward in refs. HiiJ 
contributed to this conference. 

2 Couplings of Higgs bosons to fermions 

Models with n Higgs doublets have generically a Yukawa 
Lagrangian (for the quarks) of the form: 

- £ = hf^ (/.I 4^. + h^^ 02 u'j,^ + h.c. , (2) 

where q'^, (pi, {i = 1,2) are SU(2) doublets = io'^(t>*i)i 
u'j^, d'j^ are SU(2) singlets and h"^, ft." denote 3x3 Yukawa 
matrices. To avoid flavour changing neutral couplings at 
the tree-level, it is sufficient to impose that no more than 
one Higgs doublet couples to the same right-handed field, 
as in eq. (||). 

After a rotation of the quark fields from the current 
eigenstate to the mass eigenstate basis, and an analogous 
rotation of the charged Higgs fields through a unitary 
n X n matrix U, we assume that only one of the n — 1 
charged physical Higgs bosons is light enough to lead to 
observable effects in low energy processes. The n-Higgs 
doublet model then reduces to a generalized 2HDM, with 
the following Yukawa interaction for this charged physi- 
cal Higgs boson denoted by 

jC = -^^XuL dRi + ^^^YuR^Vij dLj > H+ . 

a/2 [mw mw J 

(3) 

In (^), V is the Cabibbo-Kobayashi-Maskawa matrix 
and the symbols X and Y are defined in terms of elements 
of the matrix U (see citations in ref . El) . Notice that X and 
Y are in general complex numbers and therefore poten- 
tial sources of CP violating effects. The ordinary Type I 



and Type II 2HDMs (with n = 2), are special cases of 
this generalized class, with (X, Y) = (— cot /?, cot (3) and 
(X,Y) = (tan/3, cot/3), respectively. 

We do not attempt to list here the generic couplings 
of fermions to neutral Higgs fields. In a 2HDM of Type II, 
the neutral physical states are two CP-even states h and 
H {mil < mn) and a CP-odd state A. In this case, only 
one additional rotation matrix is needed, parametrized 
by the angle a, which is independent of the rotation an- 
gle /3 of the charged sector. This independence stops to 
be true when this model is supersymmetrized since super- 
symmetry induces a relation between tan 2a and tan 2/3. 

3 NLO corrections at the amplitude level 

The NLO corrections are calculated using the framework 
of an effective low-energy theory with five quarks, ob- 
tained by integrating out the the t-quark, the VF-boson 
and the charged Higgs boson. The relevant effective 
Hamiltonian for radiative i?-decays 

Heff = - ^ KVtb ^ a(M)0.(M) (4) 

consists precisely of the same operators Oi (/i) used in the 
SM case, weighted by the Wilson coefficients Ci{fi). The 
explicit form of the operators can be seen elsewheretl. 

Working to NLO precision means that one is resum- 
ming all the terms of the form a"(mb) ln"(mb/M), as 
well as as{mi,) (a"(m{,) ln"(m;,/M)). The symbol M 
stands for one of the heavy masses mw , fnt or mn which 
sets the order of magnitude of the matching scale 
This resummation is achieved through the following 3 
steps: 1 ) matching the full standard model theory with 
the effective theory at the scale ^iw- The Wilson coef- 
ficients are needed at the 0{as) level; 2) evolving the 
Wilson coefficients from /i = down to fi = jib, where 
fib is of the order of mf,, by solving the appropriate renor- 
malization group equations. The anomalous dimension 
matrix has to be calculated up to order a1; 3) includ- 
ing corrections to the matrix elements of the operators 
(s7|Oi(/i)|6) at the scale /i = /i;, up to order ag. 

Only step 1) gets modified when including the 
charged Higgs boson contribution to the SM one. The 
new contributions to the matching conditions Jiave been 
worked out independently by several groupsllj'ESEJ, by si- 
multaneously integrating out all heavy particles, W , t, 
and at the scale fiw- This is a reasonable approxi- 
mation provided m//± is of the same order of magnitude 
as mw or mt- 

Indeed, the lower limit on m^^i coming from LEP I, 
of 45GeV, guarantees already m^i — 0{mw)- There 
exists a higher lowet, bound from LEP II of 55GeV for 
any value of tan/3li3 for Type I and Type II models. 
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which has been recently criticized in ref. u. This criti- 
cism is based on the fact that there i^i-|no lower bound 
on niA and/or uih coming from LEP □. As it was al- 
ready mentioned, in 2HDMs (unhke in the MSSM), the 
two rotation angles of the neutral and charged Higgs sec- 
tor, a and (3, are independent parameters. Therefore, 
the pair-production process e+e^ Z* ^ hA and the 
Bjorken process e~^e~ —fZ*—^ Zh, which are sufficient 
in the MSSM to put lower limits on mh and uia sena,- 
rately, imply in this case only that mh+rriA ^ 100 GeVES. 
The other two production mechanisms possible at LEP I 
(they require larger numbers of events than LEP II can 
provide) are the decay Z h/A-f and the radiation of A 
off bb and r+r" pairs. The latter, allows for sizable rates 
only for large values of tan/? and it yields the constrapt 
tan/3 < 40 for tua = 15 GeV, in a 2HDM of Type II N. 
The former one, limits weekly_±an/3 to be in the range 
{0.2, 100} for rrih/A ~ 10 GevM. Indirect searches from 
the anomalous magnetic moment of the muon also lead 
to constraints on tan P for a light h and A. For a pseu- 
doscalar this limit is stronger than that obtained from ihe 
radiation mechanism at LEP I only for niA < 2 GeVH. 

Therefore, one of the neutral Higgs bosons can still 
be light. Charged Higgs bosons pair-produced at LEP II 
can then decay as AW^ and/or hW^ with an 

off-shell W boson. The rate is not negligible and inval- 
idates the limit of 55 GeV in Type I models and models 
of Type II with modest tan f], obtained considering only 
cs and t^Vt as possible decay modes of H^. The un- 
escapable limit of LEP I is however already large enough 
for a simultaneous integration out of H~^, W and t. 

After performing steps 1), 2), and 3), it is easy to 
obtain the quark level amplitude A{b — > 57). As the 
matrix elements {s"f\Oi\b) are proportional to the tree- 
level matrix element of the operator O7, the amplitude 
A can be written in the compact form 

A{b^s^) = ^VtbV,:D{s^\Or\b)t,ee ■ (5) 

(It should be noticed that a subset of Bremsstrahlung 
contributions was transferred to (H), as described in 
ref. Id.) For the following discussion it is useful to de- 
compose the reduced amplitude D in such a way th_at 
the dependence on the couplings X and Y (see eq. 
becomes manifest: 



XY*D, 



(6) 



In Fig. 1 the individual D quantities are shown in LO 
(dashed) and NLO (solid) order, for m//± = 100 GeV, as 
a function of /if,; all the other input parameters are taken 
at their central values, as specified in ref. El. To explain 
the situation, one can concentrate on the curves for Dxy- 
Starting from the LO curve (dashed) , the final NLO pre- 
diction is due to the change of the Wilson coefficient C7, 
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Figure 1: LO (dashed) and NLO (solid) predictions of the various 
pieces of the reduced ampUtude D for rajj± = 100 GeV (see text). 



shown by the dotted curve, and by the inclusion of the 
virtual QCD corrections to the matrix elements. This 
results into a further shift from the dotted curve to the 
solid curve. Both effects contribute with the same sign 
and with similar magnitude, as it can be seen in Fig. 1. 
The size of the NLO corrections in the term Dxy in (||) 
is 
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-40% ! 



(7) 



A similarly large correction is also obtained for Dyy- For 
the SM contribution Dsm, the situation is different: the 
corrections to the Wilson coefficient C7 and the correc- 
tions due to the virtual corrections in the matrix elements 
are smaller individually, and furthermore tend to cancel 
when combined, as shown in Fig. 1. 

The size of the corrections in D strongly depends on 
the couplings X and Y (see eq. (||)): /S.D/D is small, 
if the SM dominates, but it can reach values such as 
—50% or even worse, if the SM and the charged Higgs 
contributions are similar in size but opposite in sign. 

4 Results and Conclusions 

The branching ratio BR(i3 Xs^) can be schematically 
written as 



BR(B -> Xsi) oc |Dp -I- . 



(8) 



where the ellipses stand for Bremsstrahlung contribu- 
tions, electroweak corrections and non-perturbative ef- 
fects. As required by perturbation theory, \D\'^ in eq. 
(S) should be understood as 



1 + 2Re 



AD 



D 



LO 



(9) 



3 




m 

X 

A 
I 

IPQ 

Di 
m 



V ' ' ' ' 1 ' ' ' ' 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 
Mb= 2.4 


\\ 


4.8 - 


\ 


C7.D . . . 


- '\\ 








- 

^ V ^ 




- tan /S = 2 




1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 



100 200 300 400 500 600 
m„ 



Figure 2: BR{B Xs-y) for Y = 1, m^± = 100 GeV as a function 
of X, for Hi, = 4.8 GeV (solid), fi,, = 2.4 GeV (dahed) and fit = 9.6 
GeV (dash-dotted). Superimposed is the range of values allowed 
by the CLEO measurement. 
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i.e., the term jAD/D^'p is omitted. If Re{AD/D^^) 
is larger than 50% in magnitude and negative, the NLO 
branching ratio becomes negative, i.e. the truncation of 
the perturbative series at the NLO level is not adequate 
for the corresponding couplings X and Y. This happens 
also for modest values of X and Y, as it is illustrated 
in Fig. 2, where only real couplings are considered. The 
values X = 1 and X = — 1 in this figure, correspond 
respectively to the predictions of a the Type II and a 
Type I 2HDM with tan/5 = 1. 

Theoretical predictions for the branching ratio in 
Type II models stand, in general, on a rather solid 
ground. Fig. 3 shows the low-scale dependence of 
BR(i3 — > Xs"f) for matching scale fj,w = i^h^ i for 
?7i^± > 100 GeV. It is less than ±10% for any value 
of mfj± above the LEP I lower bound of 45 GeV. Such 
a small scale uncertainty is a generic feature of Type II 
models and remains true for values of tan/3 as small as 
0.5. In this, as in the following figures where reliable NLO 
predictions are presented, the recently calculated leading 
QED corrections^ are included in the way discussed in 
the addendum Ej of ref. EJ. They are not contained in 
the result shown in Figs. 1 and 2, which have only an 
illustrative aim. 

In Type II models, the theoretical estimate of 
BK{B Xsj) can be well above the experimental up- 
per bound of 4.5 X lO-'^ ( 95% C.L.),_reported by the 
CLEO Collaboration at this conference^, leading to con- 
straints in the (tan/?, m//±) plane. The region excluded 
by the CLEO bound, as well as by other hypotheti- 
cal experimental bounds, is given in Fig. 4. These 
contours are obtained minimizing the ratio BR(i3 
Xs7)/BR(& clvi), when varying simultaneously the 
input parameters within their errors as well as the two 
scales jUb and i^w- For tan/3 = 0.5,1,5, we exclude re- 



Figure 3: BR(B Xs7)in a Type II model with tan/3 = 2, for 
various values of /.if,. The leading QED corrections are included 
(see text). 
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Figure 4: Contour plot in (tan/3,m^±) in Type II models, ob- 
tained by using the NLO expression for BR(i? — > Xsf) and possible 
experimental upper bounds. The leading QED corrections are also 
included (see text). The allowed region is above the corresponding 
curves. 



spectively mji± < 280, 200, 170 GeV, using the present 
upper bound from CLEO. Notice that the flatness of the 
curves shown in Fig. 3 towards the higher end of mjj± , 
causes a high sensitivity jpf these bounds on all details of 
the calculation (see ref. El). These details can only alter 
the branching ratio at the 1% level, i.e. well within the 
estimated theoretical uncertainty, but they may produce 
shifts of several tens of GeV, in either direction, in the 
lower bounds quoted above. 

Also in the case of complex couplings, the results for 
BR(i3 — > Xsj) range from ill-defined, to uncertain, up to 
reliable. One particularly interesting case in which the 
perturbative expansion can be safely truncated at the 
NLO level, is identified by: Y = 1/2, X = 2exp(i(/)), and 
mjj± = 100 GeV. The corresponding branching ratios, 
shown in Fig. 5, are consistent with the CLEO measure- 
ment, even for a relatively small value of in a large 
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Figure 5: BR(i3 ^sj) as a function of <j>, where <j> parametrizes 
X = 2exp(i0), for Y = 1/2, m„± = 100 GeV. Solid, dashed and 
dash— dotted lines correspond to fii, = 4.8, 2.4, 9.6 GeV. The leading 
QED corrections are included (see text). Superimposed is the range 
of values allowed by the CLEO measurement. 



range of (f). Such a light charged Higgs can contribute to 
the decays of the t-quark, through the mode t H^b. 

The imaginary parts in the X and Y couplings in- 
duce -together with the absorptive parts of the NLO 
loop-functions- CP rate asymmetries in B ^ Xsj ■ A 
priori, these can be expected to be large. We find, how- 
ever, that choices of the couplings X and Y which render 
the branching ratio stable, induce in general small asym- 
metries, not much larger than the modest value of 1% 
obtained in the SM. 

We conclude with the most important lessons which 
can be extracted out of the calculation presented here. 
The high accuracy reached in the theoretical prediction 
of BR(B Xsj) at the NLO for the SM, is not a gen- 
eral feature. In spite of its potential sensitivity to new 
sources of chiral flavour violation, the B — Xsj decay 
may turn out to be unconstraining for many models, be- 
cause of the instability of NLO calculations. This is not 
only a temporary situation, since it is higly unlikely that 
a higher order QCD improvement is carried out. Never- 
theless, there are scenarios in which BR{B Xg^) can 
be reliably predicted at the NLO level as in 2HDMs of 
Type IL In these models, m[j± can then be safely ex- 
cluded up to values which depend on the experimentally 
maximal allowed value of BK{B Xs'y) ■ Today, we find 
m^f± ^ 165 GeV. It has to be stressed, however, that in 
the more general 2HDMs discussed in this article, 
can be much lighter and may still be detected as a decay 
product of the i-quark. 
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